Pair Breaking Caused by Magnetic Impurities in the High-Tc Superconductor 
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Conventional superconductivity is robust against the addition of impurities unless the impurities 
are magnetic in which case superconductivity is quickly suppressed. Here we present a study of the 
cuprate superconductor Bi2Sr2CaiCu2 08+i that is intentionally doped with the magnetic impurity, 
Fe. Through the use of our Tomographic Density of States (TDoS) technique, we find that while 
the superconducting gap magnitude is essentially unalTected by the inclusion of iron, the onset of 
superconductivity, Tc, and the pair- breaking rate are strongly dependent and correlated. These 
findings suggest that, in the cuprates, the pair-breaking rate is critical to the determination of Tc 
and that magnetic impurities do not disrupt the strength of pairing but rather the lifetime of the 
pairs. 



Dirty superconductors are superconductors that con- 
tain significant impurities, whether they were added in- 
tentionally or not. Conventional BCS superconductors 
are robust against normal impurities, but ruined with the 
addition of just a few magnetic impurities^ like nickel or 
manganese. This difference arises because the spin flip 
that occurs when an electron scatters off a magnetic im- 
purity violates time reversal symmetry and thus breaks 
the pair, whereas simple scattering without the spin flip 
cannot break the pair^'- . The importance of pair break- 
ing scattering is clearly evident in conventional supercon- 
ductors as seen by rapid the decrease in Tc. 

Understanding the effects of impurities on the high 
temperature superconducting cuprates is even more crit- 
ical as the cuprates are a very disordered system with 
many inherent defects. One of the best ways to study the 
effects of these defects is to intentionally add more in a 
controlled manner, through the addition of impurities^i^. 
There is still some debate whether impurities reduce Tc 
by decreasing the superconducting volume^ or increas- 
ing the pair-breaking scattering''. Therefore, cuprate 
impurity studies require a good way to measure the ef- 
fects of impurities, such as pair breaking scattering, and 
due to the material's anisotropy, to do so in a momen- 
tum resolved way. Previous studies on impurities in the 
cuprates have relied on bulk measurements^'^ or position 
sensitive spectroscopiesifiJ^^, both of which must average 
over momentum space and so provide less direct infor- 
mation about the d-wave system. 

Angle Resolved Photoemission Spectroscopy (ARPES) 
is an excellent probe to address these challenges because 
it probes the band structure directly in the momentum 
domain. It is expected that adding impurities broadens 
the band by increasing scattering and this should be di- 
rectly observable with ARPES. However, as we show be- 
low, the Fe impurities are only a small subset of the scat- 
tering events and so obtaining quantitative information 
about the Fe impurities is difficult with direct ARPES 



spectra. In the present paper we show, using our to- 
mographic density of states (TDoS) technique, that we 
can separate the pair-breaking from non pair-breaking 
scattering and extract the pair-breaking induced by Fe 
impurities, which is 1.25 meV/(% Fe). 

Here we present a study where we dope 
Bi2.iSri 9CaCu208+(5 (Bi2212) with the magnetic 
impurity^^, Fe. A series of single crystal samples of 
Bi2212 with varying concentrations of Fe were grown 
using the floating zone method^'^. Six doping levels 
were studied, with Fe concentration ranging from 0% to 
2.2%. Each sample's Tc was measured using SQUID 
magnetometry. We chose Fe impurities because they 
are known to substitute for Cui^ in the Cu02 plane. 
Consequently, the impurity potential is poorly screened 
and therefore a more significant perturbation than 
out-of-plane impurities. Increasing Fe suppresses Tc 
from 91K at 0% Fe to 67K at 2.2% Fe, a 26% decrease 
over the full range of doping (inset to FIG. The 
Tc suppression from Fe is stronger than that reported 
from either non-magnetic Zinc^"* or magnetic Ni — 
impurities, suggesting Fe has unique physics in this class 
of Cu-substituting impurities. 

Figure [T] shows a compilation of the raw nodal data 
from various doping levels. All data were taken with 
7eV photon energy and a hemispherical electron analyzer. 
The total experimental energy resolution was measured 
to be 4.5 meV using a lOK Au Fermi edge. All spectra 
were taken cold (T<30K) to minimize the effects of ther- 
mal scattering processes. One might worry that adding 
too much Fe affects the crystal structure to the point 
where a band is unrecognizable. As shown in panels al- 
a3, at all concentrations there is very clearly still a single 
band present, with very similar dispersions, widths, etc. 
Therefore, within this Fe concentration range the Fe im- 
purities are only a weak perturbation to the electronic 
structure. 

One of the standard methods to analyze ARPES data 
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the nodal spectral weight. More details about the TDoS 
method can be found in Refs j — 

The TDoS are fit to a modified Dynes formula^'' that 
includes resolution effects. Dynes' formula, Eq. ([la|. is a 
lifetime broadened BCS density of states. This formula 
was originally used in tunneling experiments on s-wave 
superconductors, where the gap is single valued over all 
of momentum space. It has been used in the cuprates in 
bulk transport and STM studies^ii^ but requires a care- 
ful integration over the d-wave gap. However, our chosen 
experimental geometry (inset to FIG. [5^2) allows us to 
treat the gap as single valued. Fitting to the adapted 
Dynes form, shown in Eq. (|lb[) . we extract the pairing 
strength. A, and the pair-breaking rate, TtdoS, from 
the TDoS spectrum. 



FIG. 1. (Color online) Panels al-a3 show raw nodal ARPES 
of Bi2212 across the measured Fe concentrations. Each spec- 
trum was taken in FY orientation at the node and at low 
temperature. Panel b the energy dependence of the MDC 
widths for all the doping levels studied. Panel c shows the 
MDCs at Ep for different dopings. Panel d shows the Im E 
extracted from the MDC at Ep for every individual sample 
studied. 



is momentum distribution curve (MDC) analysis^^. The 
MDC width is directly related to the imaginary part of 
the electronic self-energy, E"— . However, this self-energy 
includes many different types of electron scattering pro- 
cesses, and it is hard to deconvolve the individual con- 
tributions from these different processes. In particular, 
the pair-breaking scattering events are a small subset of 
the total events, so increases in this rate may be indis- 
cernible using the standard MDC analysis method. In- 
deed, as FIG. [T] panels b-d show, there is no clear trend 
in the MDC width with Fe concentration. Raw MDCs 
at Ep, in panel c, are roughly identical and the depen- 
dence of the width with binding energy, in panel b, shows 
no trend beyond the sample-sample variation. Lastly, 
panel d summarizes these findings, showing that the Im 
I](w — 0), extracted from the MDC at Ep, has no signif- 
icant trend with Fe concentration. We propose that the 
contributions from cleave-to-cleave variations in surface 
qualitjii^ overwhelm the intrinsic changes in MDC width 
from the addition of Fe impurities. 

To overcome the shortcomings of MDC analysis we 
have developed another technique to study electron inter- 
actions, the TDoS method, which is the density of states 
for a single slice through momentum space^^. (Tomogra- 
phy is the imaging of a volume via individual slices.) The 
TDoS allows us to quantitatively measure both the gap 
magnitude and the pair-breaking scattering rate. Briefly, 
to create a TDoS we isolate the coherent states of the 
band from the incoherent backgroundiS. By integrating 
the resulting spectrum across momentum, we obtain the 
coherent spectral weight. Finally, to remove effects of 
the Fermi distribution we divide this spectral weight by 
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Here /(w) is the Fermi function, i?(i?Rcs) is an energy res- 
olution term (4.5 meV FWHM Gaussian), and we leave 
the subscript TDoS on the TtdoS to distinguish it from 
a r that may be determined from EDC or MDC analysis. 

Earlier, we showed^ that the TDoS method is ro- 
bust against sample to sample variations in MDC widths 
within the same doping level. Furthermore, TtdoS is 
up to an order of magnitude smaller than the scatter- 
ing rate found via MDC analysisiS, which is consistent 
with TtdoS representing only a subset of the total scat- 
tering processes. For example, any forward scattering 
events should contribute to T," but should not signifi- 
cantly increase pair-breaking and so would not be part 
of TtdoS^^- Since the MDC is more sensitive to the to- 
tal scattering rate, it is also more sensitive to surface 
effects, such as any damage arising from cleaving, and is 
therefore a less intrinsic measure of the material. As we 
show in FIG. [21 the TDoS method can extract the pair- 
breaking dependence on Fe concentration that is simply 
masked in the MDC analysis. 

Figure [2] shows a compilation of the TDoS data for se- 
lected dopings, along with the extracted A and TtdoS 
values. Panels al-a3 show the angular dependence of the 
TDoS along with the Dynes fits for three different Fe con- 
centrations. Note how the pile-up of states increases in 
size and location with 9, signifying an increase in the gap 
size. Note also that as TtdoS increases more pairs are 
broken and more weight is present at Ep. This trend can 
be seen in the 9 = 2.5° curves (red online) moving from 
panel al to a3. One can see good quantitative agreement 
between the data and the fits with perhaps slightly de- 
graded agreement at the high angle data. We attribute 
the disagreement to a higher order effect that is not cap- 
tured by the simple nature of the functional form (which 
is only a two parameter fit). Shown in panels bl-b3 are 
the A and TtdoS values for all angles at the selected 
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Fe concentration as well as the fit to the angular depen- 
dence. A rises linearly away from the node, as a d-wave 
gap in the near nodal regime. We extrapolate this near 
nodal data to obtain the maximum gap, Amax, by fitting 
to a d-wave form: A{9) — AmaxI sm{29)\. In contrast to 
A, TtdoS is roughly constant across all angles, which is 
consistent with our previous TDoS study on BSCCO and 
is applicable over a wide doping rang o^^i^^ , but differs 
markedly from the result using MDC or EDC analysiai^. 
Due to the angular invariance, we average T^doS over all 
angles to extract Tq. 



though the TDoS technique is more quantitative than 
EDC analysis. Both of these results from Ref. 12 are 
qualitatively similar to those reported here and suggest 
that the physics involved in both Zn and Fe impurities 
is similar. Note also that while A^ax didn't change, Tq 
decreased by 26%, indicating a departure from BCS su- 
perconductivity where these quantities are directly pro- 
portional. In addition, the fact that Amax is essentially 
unchanged but both Tc and Tq change indicates that 
pair-breaking processes are critical in setting Tc in this 
material. 
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FIG. 2. (Color online) A selection of TDoS spectra for differ- 
ent Fe concentrations and the corresponding TDoS fit results. 
Panels af-a3 show TDoS curves (open circles) and the individ- 
ual TDoS fits (black lines) for selected angles from the node. 
The inset in a2 shows the k-space cuts and our definition of 
9. Panels bf-b3 show the A and TtdoS values exacted from 
the TDoS fits. The fit to A is a d-wave gap with the node 
and gap maximum as the only fitting parameters. The fit to 
FtdoS is a simple average. 

Figure|3]shows both Amax and Fq vs. Fe concentration. 
Fo rises by about 130% over the full experimental range 
while Amax remains relatively unchanged. More specif- 
ically, Fq shows a linear increase with Fe concentration, 
which we discuss in greater detail below. Despite a slight 
downward slope visually, the A^ax data show no real 
trend with added Fe; the uncertainty in the slope of a 
weighted linear fit is several times larger than the slope's 
value. Therefore, we conclude that there is no significant 
change in the superconducting gap from the addition of 
Fe impurities. A related experiment on Zn-substitutcd 
Bi2212 found, using the standard EDC method, that the 
antinodal gap size was insensitive to the addition of Zn 
impuritiesi^. In addition, Ref. 12 found a qualitative in- 
crease in the number of in-gap states with the addition 
of Zn impurities, consistent with an increased TtdoS: 
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FIG. 3. (Color online) Fq and Amax vs. Fe impurity con- 
centration. Panel (a) shows Fq, which is extracted from the 
angular average of TtdoS for a given Fe concentration. Fq 
shows a linear dependence on Fe concentration, as shown by 
the dashed line. Panel (b) shows Amax, which is extrapolated 
from the d-wave gap fit shown in figure 2. Amax shows very 
little dependence on Fe concentration, decreasing less than 
1% according to the linear fit (dashed line). The inset shows 
Tc vs Fe concentration, which changes from 91K to 67K, a 
26% change over the full Fe range. 

We propose a simple model based on ballistic transport 
that captures the essential physics of the linear increase 
of Fq with Fe concentration. The dimensionality of Fq 
suggests the simplest form would be a velocity divided 
by a length. The natural choices are the Fermi velocity 
and the mean distance between Fe impurities, respec- 
tively. Here we use the measured, nodal Fermi velocity, 
1.9 eV-A, instead of the bare velocity because the impu- 
rities are a perturbation on a system that already has 
self-energy effects. In this simple model, we estimate the 
mean distance as the Cu-Cu spacing, 3.82 A, divided by 
the concentration, x, of Fe scattering sites. And we add a 
dimensionless parameter, a, that represents the 2D scat- 
tering cross-section for these events. Note that in two 
dimensions, the scattering cross-section has dimensions 
of length instead of the normal dimensions of area for 
a 3D scattering process. Since we have divided out the 
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length dimension using "a" , a is dimensionless and has 
units of Cu-Cu spacing. This model would be sufficient 
if Fe impurities were the only things causing scattering. 
However, Bi2212 has a number of native defects includ- 
ing Cu vacancies and several types of oxygen vacancies 
and dislocations^^"— . Each of these scattering species 
should contribute to the scattering rate in the same way 
as outlined for Fe impurities, and so we can add a sum 
of terms to the formula for Fq, as shown in Eq. [2al This 
sum includes all native defects in the material, including 
any not listed above. However, it is beyond the scope 
of this paper to determine both the concentrations and 
cross-sections for all these defects. And since they are 
not relevant to extracting properties of the Fe impuri- 
ties, we group all of the native defects together in a term 
called F^^*'™, shown in Eq.[2bl This model should hold 
as long as the impurities are well separated and therefore 
not interacting with each other. 
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Figure 2] shows the linear fit to our data in the con- 
text of this model. The linear increase in Fq we at- 
tribute to Fe impurity scattering while the offset is due 
to the combined effects of all native defects in the ma- 
terial. From this fit we extract a = 0.25, which means 
the pair-breaking scattering cross-section of a single Fe 
impurity is approximately 25% of the Cu-Cu spacing, or 
0.96 A. This is a reasonable value of the cross-section 
for a single site defect, and suggests that Fe impurities 
are relatively weak pair-breakers. This result is also con- 
sistent with STM work on Fe-substituted BSCCO that 
claims Fe is a weaker scatterer than NiiSi2^, Znii, and 
Cu vacancies^,—. To our knowledge, this is the first 
time the pair-breaking cross-section has been measured 
for any impurity in the cuprates. Furthermore, this gen- 
eral experimental procedure would work on any impurity, 
as long as there is sufficient control to increase the con- 
centration of a known impurity while keeping all others 
constant. 

Using the TDoS technique we have measured the ef- 
fects of Fe impurities on the superconducting gap and 
pair-breaking scattering rate in the high Tc supercon- 
ductor Bi2212. The main effect of Fe impurities is to 
increase the pair-breaking scattering rate, while leaving 
the gap parameter unchanged. This indicates that mag- 
netic impurities do not affect the pairing strength, just 



the pair lifetimes. This result also confirms that cuprate 
superconductivity is sensitive to pair-breaking scatter- 
ing. The correlation between Fq and Tc is suggestive 
and motivates further studies at temperatures near Tc, 
as well as studies of other impurities, such as Zn and Ni, 
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FIG. 4. (Color online) Segregating the effects of Fe impuri- 
ties from native defects using our simple model of impurity 
scattering, pjj'''*'™ is the pair breaking rate for all the com- 
bined native defects not including any Fe impurities, while a 
is a dimensionless parameter that represents the pair-breaking 
scattering cross-section (in units of Cu-Cu spacing) for Fe im- 
purities. 



in hopes of discovering any universal behavior of cuprate 
impurities. 
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